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Abstract 

The large-scale dynamics of matter is inferred from the observed peculiar velocities 
of galaxies via the POTENT procedure. The smoothed fields of velocity and 
mass-density fluctuations are recovered from the current data of ~ 3, 000 galaxies. 
The cosmological density parameter O can then be constrained in three ways: (a) 
by comparing the density fields of mass and galaxies, (b) by using the velocity field 
in voids, and (c) by investigating quasilinear deviations from Gaussian fluctuations. 
The results indicate a high value of fi ~ 1; values in the range 0.1-0.3 are rejected 
with high confidence. 



1. Introduction 

The raw data are estimated distances V{ and redshifts Z{ for a set of objects in directions 
fi [34] The distances are estimated by the Tully-Fisher method for spirals, and by the 
analogous D n — cr method for ellipticals and SO's. The corresponding radial peculiar velocities 
are m = czi — (H is set to unity; distances are measured in kms -1 ). Given these 
sparsely-sampled radial velocities, POTENT first computes a smoothed radial- velocity field, 
u(r), in a spherical grid using a tensor window function [8,10]. In the current analysis we use 
a Gaussian window with radius ~ 1200 kms -1 . Weighting inversely by the local density near 
each object mimics equal- volume averaging which minimizes the bias due to the nonuniform 
sampling, and weighting inversely by the distance variance, af , reduces the random effects of 
the measurement errors. 

The velocity field is recovered under the assumption of potential flow: vix ) = — V$(as) [3]. 
According to linear gravitational instability theory (GI) any vorticity mode decays in time 
as the universe expands, and based on Kelvin's circulation theorem the flow remains 
vorticity-free in the mildely-nonlinear regime as long as the flow is laminar. The velocity 
potential can therefore be calculated by integrating the radial velocity along radial rays, 
= — J Q r u(r' , 6, cf))dr' . Differentiating $ in the transverse directions then recovers the 
two missing velocity components. 



The underlying mass-density fluctuation field, 8(x) } is computed by the approximation [28] 



8 c (x) = \\I - /(ft) _1 dw/da?|| - 1 , (1) 

where the bars denote the Jacobian determinant, / is the unit matrix, and /(O) = D / D ~ 
O - 6 with D(t) the linear growth factor [30]. Eq. 2 is the solution to the continuity equation 
under the Zel'dovich assumption that particle displacements evolve in a universal rate [35]. 
This expression still involves only the first partial derivatives in Eulerian space, and it reduces 
to the familiar 8 = — /( Q )~ ' V • v in the linear regime. It has been found to approximate the 
true density in N-body simulations with an rms error less than 0.1 over the range —0.8 < 
8 < 4.5 [26,23]. 

The results suffer from random and systematic errors due to the distance errors and the 
nonuniform sampling. The ways we handle these errors are discussed in detail in [8] and 
in [10], where special attention has been paid to correcting for inhomogeneous Malmquist 
bias. 

Figures 1 and 2 show preliminary maps of the recovered density field, now extending to 
beyond 8000 kms -1 in some regions [10]. The extended data includes the sample of 493 
objects [6,22] which was the basis for the original POTENT90 analysis [9], and several much 
bigger spiral samples [7,15,16,24,25,33] - a total of more than 3000 galaxies that we have 
self-consistently calibrated and carefully put together [12,34]. 
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Figure 1: Cosmography: the fluctuation fields of velocity and mass-density in the 
Supergalactic plane as recovered by POTENT from the velocities of ~ 3000 galaxies 
with 1200 kms -1 smoothing. The vectors shown are projections of the 3D velocity 
field in the CMB frame. Contour spacing is 0.2 in 8 } with the heavy contour marking 
8 = and dashed contours denoting negative fluctuations. The Local Group is at 
the origin. The GA is on the left, PP on the right, and Coma is at the top. 

The density peak of the Great Attractor (GA) under 1200 kms -1 smoothing is of 8 = 1.2±0.3 
near the Galactic plane (X ~ —4000, Y ~ 0). The Perseus-Pisces (PP) peak is at the level of 
8 = 0.7 ±0.3 and it extends towards Aquarius and Cetus near the south Galactic pole (where 
the "Southern Wall" is seen in redshift surveys). The Coma Great Wall shows up towards 
the north galactic pole. Two great voids extend from bottom left to top right. 
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Figure 2: Cosmological landscape: mass density fluctuations in the Supergalactic 
plane as recovered by POTENT from the velocities of ~ 3000 galaxies with 
1000 kms -1 smoothing. The height of the surface is proportional to 8. The map 
extends to a distance of 8000 kms -1 about the Local Group. The GA is on the left, 
PP on the right, and Coma at the top. 

Is there a back flow behind the GA? Is PP moving toward the LG or away from it? These 
effects are detected by the current POTENT analysis still only at the 1.5-sigma level. Recall 
that there is an important free parameter in the velocity field: the distances are determined 
only relative to each other and the mean Hubble flow is not known. It has been determined 
here by minimizing residuals within the volume sampled, but a Hubble-like peculiar velocity 
error of order 10% is not out of the question. The data thus allow a solution with no GA 
back flow and with PP moving away from the LG. 

A simple and robust statistic is the bulk velocity within a sphere of radius R centered on the 
LG. It is computed by weighted vector averaging of the smoothed peculiar velocities at the grid 
points. The results suffer from a systematic uncertainty which arises from the non-uniform 
sampling. We show in Figure 3 two different results, obtained by applying different relative 
weights to the data; one minimizing the sampling-gradient bias by equal-volume weighting, 
and the other minimizing the random errors by weighting with the inverse of the distance 
variances. The difference between the two results reflects the systematic uncertainty. The 
bulk velocity in a top-hat sphere of radius 6000 kms -1 is thus in the range 270 — 360 kms -1 . 
The additional 1-sigma random error in \V\ due to distance errors is typically 15%, as derived 
from Monte-Carlo noise simulations. This error does not include the random error due to 
cosmic scatter - the fact that only one sphere has been sampled. 

Three arguments support our belief that the peculiar velocities are real, and consistent with 
being generated by GI: (a) The velocity fields traced separately by spirals and by ellipticals, 
using different distance indicators, are consistent with each other [21] (b) the COBE detection 
of 8T/T ~ 10 -5 and the bulk velocity of ~ 350 kms -1 measured in the local neighborhood in 
a top-hat sphere of radius 6000 kms -1 are consistent with each other [5], and (c) the galaxy 
density from redshift catalogs and the mass density recovered from observed velocities are 
consistent with each other (§2; [9]). We thus assume that the velocities are real and generated 
by the gravity of the mass-density fluctuations, and try to determine O. 
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Figure 3: The bulk velocity in a top-hat sphere of radius R e ff about the LG. 
Shown are \V\ (filled), V x (triangles), V y (squares), and V z (hexagons). The two 
results shown reflect the systematic uncertainty. The 1-sigma random error due to 
distance errors is ~ 15%. 

2. Galaxies Versus Mass: Q and Biasing 

The POTENT density 8 P} determined assuming = 1, relates to the true 8 in the linear 
regime by 8 P oc f(Q)8. On the other hand, if we parameterize the relation between the galaxy 
density fluctuation field, 8 G} and the mass density fluctuations by a universal "biasing" factor, 
S g = bS } then we expect a relation of the sort 8 P = [f(Q)/b]8 G . Given the uncertainties in the 
two datasets, we can ask whether the POTENT data are consistent with being a noisy version 
of the galaxy data, and obtain the best-fit value for f(Q)/b with associated confidence limits. 

The degeneracy of O and b is broken in nonlinear regions, where 8(v) is no longer simply 
proportional to /(fi) _1 . The quasilinear analysis based on Eq. 1 allows a first attempt at 
determining O and b separately. 
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Figure 4: Galaxies versus mass density fields in the Supergalactic plane. Contour 
spacing is 0.2 in 8, with the heavy contour marking 8 = 0. The Local Group is at 
the origin. GA is on the left, PP on the right, and Coma at the top. 



Figure 4 compares the density maps in the Supergalactic plane of the POTENT mass density 
and the galaxy density from a redshift survey of IRAS galaxies, flux limited at 1.2Jy [13], both 



Gaussian smoothed at 1200 kms -1 , assuming 6=1. Given the errors, the similarity between 
the maps is strong. Both feature the general GA phenomenon as a ramp which peaks beyond 
the Hydra-Centaurus clusters (X ~ —4000) and falls off gradually toward Virgo (X ~ —300, 
Y ~ 1300) and toward Pavo across the Galactic plane. The fields also agree well in the PP 
and Cetus superclusters, in the Coma Great Wall region (X ~ 0, Y ~ 8000), and in the great 
void in between. 

These new data have not been subject to a quantitative comparison yet. So far, we have 
applied an elaborate statistical analysis [9] only to the earlier data of POTENT90 and IRAS 
1.9Jy [31], Gaussian smoothed at 1200 kms -1 . Noise considerations limited that analysis to 
a volume ~ (5300 kms -1 ) 3 containing ~ 12 independent density samples. Monte-Carlo noise 
simulations showed that the data are consistent with the hypotheses of linear biasing and GI. 
The Monte-Carlo simulations were then used to estimate the random errors, to correct for 
systematic errors in POTENT, and to constrain the parameters via a likelihood analysis. Our 
robust result is fi ' 6 /6 7 = 1.28lo';5g a ^ 95% confidence. Small nonlinear effects allow weaker, 
separate constraints on O and on 6 7 : if ft = 1 then 6 7 = Q.7±™, and if 6, > 0.5 then O > 0.46. 
Inhomogeneous Malmquist bias could decrease these estimates of O, but our 95% confidence 
limit for 6 7 > 0.5 could be reduced at most to O > 0.3. The range of uncertainty is big 
because the data were limited and because we were extremely careful in our error analysis. 
The preliminary results from a similar comparison with the galaxy distribution in the optical 
catalog [17,18] indicate a similar correlation between light and mass, with /(0)/6 o ~ 0.7±0.2, 
in agreement with the ratio of 6 /6 7 ~ 1.5 obtained by direct comparison, and in agreement 
with common theoretical ideas. Quantitative comparisons using the extended velocity data 
and the new IRAS and optical surveys will hopefully be able to provide clearer answers to the 
fundamental questions concerning the validity of GI, the biasing scheme, and the value of O 
with a reduced range of uncertainty. 

3. Q from Voids 

A diverging flow in an extended low-density region can provide a robust dynamical lower 
bound on O, based on the fact that large outflows are not expected in a low-fi universe [11]. 
The velocities are assumed to be induced by gravity from small initial fluctuations, but no 
assumptions need to be made regarding their exact Gaussian nature, galaxy biasing, or A. 
The derivatives of a diverging velocity field infer a nonlinear approximation to the mass 
density, 8 C (Eq. 1), which is an overestimate when the true value of O is assumed []. 
Analogously to the behavior of the linear approximation, So = — Q -(U: V • v, the 8 C inferred 
from a given velocity field is a monotonically increasing function of O, and it may become 
smaller than —1 for O values that are too small below the true value. Then, since 8 > — 1 
because mass is never negative, the allowed values for O are bounded from below. 
Given the observed radial peculiar velocities of galaxies, one can use the POTENT procedure 
to recover the three-dimensional velocity field, Gaussian smoothed at ~ 1200 kms -1 . One 
can then derive the inferred mass-density field and the associated error field for different 
values of O. Focusing on the deepest density wells, the assumed O should be lowered until 
8 C becomes significantly smaller than —1, where the significance is estimated based on the 
estimated errors. Such low values of O could be ruled out. 

The most promising "test case" provided by the current data seems to be a broad diverging 
region centered near the supergalactic plane at the vicinity of (X, Y) = (—2500,-4000) in 
kms -1 supergalactic coordinates. This region can be roughly identified with the "Sculpter 
void" of galaxies, clearly seen in the SSRS redshift catalog next to the "Southern Wall". 
As can be seen in Figure 5, values of O 1 are perfectly consistent with the data, but 8 C 




X (1000 km/sec) X (1000 km/sec) 

Figure 5: Maps of the density-fluctuation field inferred from the observed velocity 
derivatives, 8 C} in part of the supergalactic plane, for two different values assumed 
for 0. The void of interest is confined by the Pavo part of the Great Attractor on 
the left and the Aquarius extension of the Perseus-Pisces supercluster on the right. 
The Local Group is marked by a '+'. Contour spacing is 0.5, with the mean, 8 C = 0, 
marked by a heavy line, 8 C > solid, and 8 C < dotted. The heavy-dashed contours 
mark the downward deviation of 8 C from —1 in units of the standard deviation <r, 
starting from zero (which coincides with 8 C = —1), and decreasing with spacing 
—0.5cr. The value O = 0.2 is ruled out at the 2.9-sigma level. 

becomes smaller than —1 in this void already for O = 0.6. The confidence by which 8 C < — 1 
at a given point is expressed in terms of the random error a there. The values O = 0.4, 0.3, 
and 0.2 are ruled out at the 1.6-, 2.4-, and 2.9-sigma levels respectively. 

The result reported here is just a preliminary application of the method. The systematic errors 
have been partially corrected for in the POTENT procedure, but a more specific investigation 
of the biases affecting the smoothed velocity field in the deepest density wells is required for 
more reliable results [14]. 

For the method to be effective we need to find a void that is (a) bigger than the correlation 
length for its vicinity to represent the universal O, (b) deep enough for the lower bound to 
be strong, (c) nearby enough for the distance errors to be small, and (d) properly sampled 
to trace the velocity field in its vicinity. 

4. from the Probability Distribution Function 

A generalization of the method using voids to determine O makes use of the whole one-point 
probability distribution function (pdf) of the smoothed POTENT density field. Assuming that 
the initial fluctuations are a random Gaussian field, the PDF dvelops a characteristic skewness 
due to nonlinear effects, which is noticieable early in the quasilinear regime [23]. The skewness 
of 8, in second-order perturbation theory, is approximated by (8 3 )/(8 2 ) 2 = (34/7 — 3 — n), with 
n the effective power index near the smoothing scale [24]. This result is almost independent 
of O. Since in linear theory 8 = — /(fi) -1 V-u, one can expect the corresponding skewness of 
V • v to strongly depend on O. Indeed, our second-order calculation yields [25] 

<> 3 = ((V • t0 3 )/((V • v) 2 ) 2 = -/(0)" 1 (26/7 - 3 - n) . (2) 

Using N-body simulations and 1200 kms -1 smoothing we indeed find S3 = 1.8 ±0.6 for 0=1 
and S3 = 4.1 ± 1.1 for O = 0.3. The quoted error is the standard deviation associated with 



computing S3 within one single sphere of radius 5000 kms 1 in a CDM universe (Hq = 75, 
6=1). 

The value of S3 in the current POTENT velocity field within 5000 kms -1 is 0.6 ± 0.5, here 
the error representing distance-measurement errors. The two kind of errors should roughly 
add in quadrature. Hence, O = 0.3 is rejected at the 2.5-sigma level. This bound is similar 
in strength to the bound obtained using voids, but it relies on the assumption of Gaussian 
initial fluctuations and it is somewhat sensitive to the exact shape of the power spectrum 
(mostly through the error estimate). 

The method just described still makes use of only a fraction of the data - the PDF - which 
does not fully specify the field. In fact, the PDF tends to develop a general lognormal 
shape in the course of quasilinear evolution even in certain cases of non-Gaussian initial 
fluctuations [26]. A more powerful bound on O can be obtained if one makes use of the 
detailed POTENT velocity field to recover the initial PDF (jpdf), and then use this I PDF 
to constrain O. Nusser and Dekel [27,28] have developed a method for tracing quasilinear 
fluctuations back in time and determining the IPDF using an Eulerian representation of the 
Zel'dovich approximation. Then, the key point for determining O is that the IPDF recovered 
from the present velocities is very sensitive to the value of O assumed in the recovery process 
(similar in principle but opposite in sign to the O dependence of the present PDF of V • v). 
If we are willing to assume a specific shape for the IPDF based on other considerations, most 
naturally a Gaussian shape, then we can constrain O by comparing the recovered IPDF with 
the assumed one. 
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Figure 6: The density IPDF recovered from the velocity field provided by 
POTENT90 from observed velocities (solid) compared with Gaussian (short dash) 
and with the IPDF recovered from the velocity field of Gaussian CDM simulations 
(triangles). The assumed O is 1 or 0.3. The simulations are of fi = ^ accordingly. 

The velocity field out of POTENT90 within a conservatively selected volume has been fed 
into the IPDF recovery procedure with O either 1 or 0.3. The recovered lPDF } s are shown in 
Figure 6. The error bars attached to the POTENT IPDF are the Monte-Carlo measurement 
errors, and the error bars attached to the model IPDF estimate the error due to the limited 
volume sampled. The total error is roughly a sum in quadrature of the two errors. The 
immediate impression is that the IPDF recovered with O = 1 is marginaly consistent with 
Gaussian while the one recovered with O = 0.3 significantly deviates from a Gaussian. Trying 
to reject these models we use several different statistics which characterize the IPDF. The 
standard deviation of each statistic due to measurement and volume errors was evaluated 
using Monte-Carlo simulations. For example, bin by bin in the IPDF, the largest deviation 



for fi = 1 is ~2-sigma while for fi = 0.3 there are deviations larger than 4-sigma near 
x = (8 — /i)/cr ±1.5. The standard moments of skewness and kurtosis are poorly determined 
because they are tail-dominated, but the replacements and (\x\) do very well: The 

hypothesis fi = 0.3 is rejected at the 6.2-sigma level by the former and at 5.3-sigma by the 
latter, while fi = 1 is "rejected" at below the 1.6-sigma level by both. Thus, fi > 0.3 is 
strongly ruled out! 

The bounds on fi discussed in this section depend on the assumption of Gaussian ipdf. But 
in fact, the IPDF does not have to be assumed: it can be determined from an observed density 
field independently of fi [27]. We have recovered the IPDF from the 1.2Jy IRAS density field 
within 80 h -1 Mpc of us, assuming 6 7 = 1. The method is not that sensitive to the exact 
value of bj in the range 0.5 — 1.5. The reconstructed IPDF is consistent with Gaussian at the 
~l-sigma level even without taking into account the measurement errors in the density field 
itself [29]. 

5. Conclusions 

The maps of velocity and mass-density fields obtained from extended data of ~ 3000 galaxies 
with 1200 kms -1 smoothing show new dynamical features such as the Perseus-Pisces ramp 
extending to Aquarius and Cetus, the Coma Great Wall and the great voids in between. The 
bulk velocity within a top-hat sphere of radius 6000 kms -1 is in the range 270 — 360 kms -1 . 

Based on the recovered dynamical fields one can put the following constraints on fi: 

• By comparing POTENT90 mass density with 1.9Jy IRAS galaxy density: O - 6 /6 7 = l^l^g 
at 95% confidence. The extended data indicate a value near unity with tighter constraints. 

• A comparison with optical galaxy density indicates fi ' 6 /& o ~ 0.7. These results are 
consistent with fi = 1, b T = 1 and b a = 1.5 at ~ 1200 kms -1 ; a much lower value for fi would 
require unrealistic antibiasing for IRAS galaxies. 

• The velocity around the great local void requires fi > 0.3 at the 2. 4-sigma level, independent 
of biasing. 

• If we adopt the hypothesis that the initial fluctuations were Gaussian, supported by our 
finding that the density of IRAS galaxies unambiguously indicates a Gaussian IPDF, we find 
by the skewness of the distribution of V • v that fi > 0.3 at the 2.5-sigma level. 

• Under the same assumption of Gaussian IPDF, the IPDF recovered by POTENT90 from the 
observed velocities says that fi > 0.3 at the 6-sigma level 

The velocity data thus indicates a high value for fi, near unity, with fi < 0.3 strongly ruled 
out. The uncertainty in this general conclusion is mostly due to inhomogeneous Malmquist 
bias in the velocity data. Tests of our correction method indicate that this bias cannot affect 
the resultant fi by more than 20%. 
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